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Abstract 

Using the operator product expansion we show that the 0(1/ p 2 ) correction to the pertur- 
bative expressions for the gluon propagator and the strong coupling constant resulting from 
lattice simulations in the Landau gauge are due to a non- vanishing vacuum expectation value 
of the operator A^A^. This is done using the recently published Wilson coefficients of the 
identity operator computed to third order, and the subdominant Wilson coefficient computed 
in this paper to the leading logarithm. As a test of the applicability of OPE we compare 
the (A^A^) estimated from the gluon propagator and the one from the coupling constant 
in the flavourless case. Both agree within the statistical uncertainty: J (A^A^) ~ 1.64(15) 
GeV. Simultaneously we fit Aj^g- = 233(28) MeV in perfect agreement with previous lattice 
estimates. When the leading coefficients are only expanded to two loops, the two estimates 
of the condensate differ drastically. As a consequence we insist that OPE can be applied in 
predicting physical quantities only if the Wilson coefficients are computed to a high enough 
perturbative order. 

P.A.C.S.: 12.38.Aw; 12.38.Gc; 12.38.Cy; 11.15.H 



1 Introduction 

When computing in a fixed gauge an operator product, the operator product expansion 
(OPE) contains in general contributions from local gauge-dependent operators, even though 
they should not emerge in the gauge-invariant sector. For example in ref. Q, a detailed analysis 
clearly shows that operators such as A 2 = A^A^ contribute to QCD propagators' OPE through 
a non-zero expectation value in a non-gauge-invariant "vacuum". A 2 is the unique dimension 
two operator allowed to have a vacuum expectation value (v.e.v.) and is thus the dominant non 
perturbative contributor, leading to ~ 1/p 2 corrections to the perturbative result. 

These expected ~ 1/p 2 have at first sight nothing to do with the possible presence of 1/p 2 
terms in gauge invariant quantities such as Wilson loops 0: since no local gauge invariant 
gluonic operator of dimension less than 4 exists it is expected from OPE that the dominant 
power correction should be cx 1/p 4 , originating from the local and gauge-invariant G^ U G^ W . Of 
course the operator A 2 in the Landau gauge can be viewed, by simply averaging it over the gauge 
orbit, as a gauge invariant non local operator. But then, dealing with non-local operators, we 
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loose the standard OPE power counting rule relating the power behaviour of a Wilson coefficient 
to the dimension of the corresponding operator: there is no reason for this non-local operator 
to yield 1/p 2 contributions in a gauge invariant observable. It has been strongly stressed in 
ref. [|| that, working in the Landau gauge, the A 2 operator plays a special role since, imposing 
the Landau gauge condition is equivalent to asserting that A 2 is at an extremum or a saddle 
point on its gauge orbit. Practically, on a lattice, one fixes the Landau gauge by searching 
for a minimum of A 2 on the orbit. We are not able to elaborate further on the issue of what 
relation might exist between the expected (A 2 ) condensate in the landau gauge and the possible 
unexpected 1/p 2 terms in gauge invariant quantities [§]. But we are in a position to put the first 
step of this possible route on a firm ground: to provide a strong evidence that there is indeed 
an (A 2 ) condensate in the Landau gauge and that it is not small. 

To that aim we will use heavily OPE. We need to be sure that OPE really works in this 
situation and have to invent some way of verifying this point. A success of this check would 
achieve several goals. First it would give a strong support to the conjecture that OPE is really 
working in this situation, i.e. that we do not encounter a strange situation where OPE would 
have failed like the one discussed in the preceding paragraph about (G^G^). Second it would 
confirm that we go far enough in the perturbative expansion (the expansion in l/ln(p 2 )) to be 
able to say something sensible about the power expansion (in 1/p 2 ). Third it would confirm 
that we really are measuring {A 2 ). Such checks have of course many consequences which will be 
further discussed in the conclusions. 

From a practical (numerical) point of view, 1/p 2 terms provide a specially convenient way 
to test OPE since 

- they remain visible at much larger energies than the 1/p 4 ones which would result from 
the gauge invariant G^G^ 

- as already mentioned their OPE analysis is rather simple and unambiguous because A 2 is 
the only dimension-two operator to contribute. 

A recent study of a s MOM (p), the Landau gauge coupling constant^, regularized on a lattice 
showed unequivocally the presence of 1/p 2 power corrections still visible at energies ~ 10 GeV 
for which OPE contributions of the gluon condensate, (A 2 ), were natural candidates In this 
term all the non-perturbative input is contained in {A 2 ) while the OPE Wilson coefficients can 
be computed in perturbation. In view of this, we proposed in a previous work || a procedure to 
test OPE based on the determination, and further comparison, of the two estimates of the gluon 
condensate, {A 2 ), obtained from both gluon two- and three-point Green functions by means of 
a simultaneous matching of the lattice data to the OPE formulas derived by following standard 
Shifman-Vainshtein-Zakharov (SVZ) techniques ||. Thus, our OPE matchings of lattice data 
provides two independent estimates of the renormalized A 2 condensate. The adequate definition 
of renormalized condensates and their "universality" when studying different Green functions 
was discussed in ref. [6f in connection with the choice of truncation orders for perturbative and 
OPE series (see also [H],[Tl]]). In this preliminary work we described the theoretical framework 
for this testing procedure and we performed a first analysis of previous lattice data |7[ || but the 
perturbative /^-function was known at that time only up to two loops and our use of OPE was 
limited to the sole computation of the Wilson coefficients of A 2 at tree-level 

After this work was completed a computation of the third coefficient of the MOM beta 
function, /?2 has been published in ref. Q. The authors of this last work conclude that their 

2 a s MOM (j>) stands for the QCD running coupling constant non-perturbatively renormalized in a kinematically 
asymmetric point by following the momentum subtraction prescriptions (MOM) 
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computed 02 and our "prediction" of this coefficient based on OPE consistency [|| reasonably 
agree with each other. Thanks to the new information concerning the (3- function and to the 
high accuracy of our lattice results we are now in a particularly favourable situation to address 
further the questions we have mentionned above. This is the task we shall attack in the present 
paper, presenting a consistent calculation in the MOM scheme (a symmetric kinematics chosen 
for the vertex) with the Wilson coefficients of the identity operator computed at three loops 
Hg, flU and the ones of A 2 computed to the leading logarithm in section |2.2| . In particular 
we will compare the check of the "universality" of the condensates when expanding the leading 
perturbative coefficients to three loops and when one uses only the two loop order. 

The theoretical setting of our use of OPE is described in Section ||: the tree-level computa- 
tion, presented previously in ref. ||, is only sketched and most attention is paid to the obtention 
of the one-loop anomalous dimension of the Wilson coefficients. The fitting strategy is explained 
and the matching test performed in ^. Finally, we discuss and conclude in Section ^. 

2 OPE for the gluon propagator and a s (p) 

In the present section we shall expand the three-point Green function, and hence as(p), as 
well as the gluon propagator, in the OPE approach up to the the 1/p 2 order. Both gluonic two- 
and three-point Green functions are renormalized according to the MOM scheme. Let us start 
with a reminder of the computation of the tree-level Wilson coefficients ||. 

2.1 Tree-level Wilson coefficients 

In the pure Yangs-Mills QCD, without quarks, OPE yields 



+ 
+ 



(1) 



(*)2f,(pi,P2,PJ») 1 + (di)^(pi,P2,Ps) : A$(0) : 

+ (d 2 ) a ^X(pi,P2,P3) :<,(0)4'(0): + ••• ; 

(2) 

where only normal products of local gluon field operators occur and A (A) stands for the gluon 
field in configuration (momentum) space, a, b being colour indices and [a, v Lorentz ones. The 
notation T() simply refers to the standard T* product in momentum space. The normal product 
of Eqs. (|ll§) should be defined in reference to the perturbative vacuum ||. 

Only terms in Eqs. (HUH) containing an even number of local gluon fields give a non-null . 
v.e.v. because of Lorentz invariance and of the gauge condition ||. The coefficients Co and do 
are the purely perturbative Green functions. Assuming the Wilson factorization of soft and hard 
gluon contributions, the relevant Wilson coefficients C2,d,2 can be obtained, in perturbation, by 
computing the diagrams in Fig. [I] which represent the matrix elements of operators in the l.h.s. 
of Eqs. (||J2|) between soft gluons, indicated by crosses. Using also the tree level expression for 
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Figure 1: Four and five gluons tree-level diagrams contributing (with all their possible permutations) to 
the Wilson coefficients of the gluon propagator and the three-gluon vertex. Crosses mark the gluon legs 
due to the external soft gluons. 



the matrix element, between the same two soft gluons, of the local operators in the r.h.s. of 
Eqs. (|l|j2|) we obtain C2 and d,2 by matching both sides. 

Thus, in the appropriate Euclidean metrics for matching to lattice non-perturbative results, 
we can write: 



k 2 G^(k 2 ) = z MOM (k 2 ) = zl™(k 2 ) +i^iy p. 

k e G^ 3 \k 2 , k 2 , k 2 ) = fc 6 Gp2-t(& 2 ) k 2 , k 2 ) + ^^y 2 _ -g I (3) 
where the scalar form factors , G^ are defined as follows from the Green functions 



G^(k 2 ,k 2 ,k 2 



N*-l 3 



abc 



(A a J-p) A b u {p)) 



-(Al{ Pl )A b MA c p (pz)) 



18/fc 2 N c (N?-l)' ^ LJ pK 



treeN/Ji^2M3 



+ 



(pi ~ P2) P (P2 ~ PzYiPZ ~ PlY 

2k 2 



(4) 



For the kinematic configuration p\ = p 2 = p 2 = k 2 the three-gluon tree-level tensor is defined as 



(T 



tree \ 



lfllfJ,2fJ-3 



Vi/4 (Pi ~ P^t^'a + cycl.perm.J ( 5^. 

i=l,3 ^ 



Pi 



(5) 



In Eqs. @-@ we have dealt with bare quantities, depending only on the cut-off a -1 and 
on the momentum k. We have omitted to explicitate the dependence on the cut-off in order to 
simplify the notations. Using Eqs. (||) these Green functions can be conveniently renormalized 



A 



by MOM prescriptions: the renormalized two-point Green function is taken equal to 1/k 2 for 
k 2 = V 2 , 



h 2 r m (h 2 2 ,_ k 2 G^(k 2 ) fk 2 \ (k 2 \ (A*)^ 1 
The Co Wilson coefficient can be written as 

* 2 , A ^o°o^ 2 ) . „ , « «(^) 



co( -^,a(/x)J = zMOM( ^ 2) = c (l,«(/i))^MOM^2) ' ( 7 ) 



MOM 



and verifies consequently the perturbative evolution equations of Z 

dlnco(jJ,a(/x)) / a ( A .) fa(k)\ 2 _ /a(*0V 



70^r +71 TT +T5 TT + • • ■ ( 8 ) 



dlnfc 2 I ,u 4vr ' V 4tt / ' V 4?r 

where 71 and 72 depend on the perturbative scheme in which the strong coupling constant a(k) is 
defined. The boundary condition to solve Eq. (||) comes from the nonperturbative normalization 
of k 2 G^\k 2 ,fJL 2 ) to 1 at k 2 = (i 2 , and it results that c (l,a(/x)) = 1 + 0(1/^ 2 ). 

Let us remind that in the MOM prescription, the three-point Green function is renormalized 
by G^ (k 2 , /I 2 ) = G^(k 2 ,k 2 ,k 2 )(Z MOM (fj,))- 3 / 2 , and the MOM coupling constant follows from 

„ fu2\ _ G^jk 2 , k 2 , k 2 ) / MOM/.2^3/2 
9R{k ] ~ ( G (2)( fc 2 )) 3 ^ ^ )) 

= fc 6 4 3) (fc 2 ,/, 2 )(fe 2 4 2) (fc 2 ,/, 2 )) -3/2 . (9) 

Analogously to Eq. (g) we define the renormalized three point Green function 

fc 6 <4% 2 , M 2 ) = do + <*2 4^T) ^2 • (10) 

where the do Wilson coefficient verifies the perturbative evolution equations of k s G$ (k 2 , fj, 2 ) 
and the boundary condition do(l, a(fj,)) = git(fJ, 2 ) + 0(l//x 2 ) is immediate from Eqs. (|g|JTo|). 

In the MOM scheme, the gluon condensate, {A 2 )r^, is renormalized at fj? by a standard 
condition, through division by a renormalization constant Z42 . 

The C2 and d% Wilson coefficients at tree level are B 



c 2 (1, «(//)) = 3 g 2 , 

da(l,a0/)) = 9 5 3 • (H) 

Since the three point Green function naturally defines the MOM scheme coupling constant 
(see below @), we will perform all the coming calculations in the symmetric MOM scheme 
where 

70 = 13/2, 71 = -16.9, 72 ~ 1332.3 . (12) 



In Eq. @ a(k) is of course taken to be the purely perturbative running coupling constant, 
5r port(^ 2 )/(^ 7r )' obtained by integrating the beta function in the MOM scheme, 

^ a(k) = P(a(k)) = - (|V(fc) + A Q 3 (fc) + _|L_ J {k) (13) 
where j|] 



A) = 11, A = 51, 



3088. 



(14) 



2.2 Wilson coefficient at leading logs 

The purpose is now to compute to leading logarithms the subleading Wilson coefficients 
in Eqs. (|6,10). To this goal, following Jl4| it will be useful to consider the following matrix 
element, 



{g a T \A r p {k)A%{-k)\g b u ) R ^ 



'pa 



kp k a 



k A 



4(N 2 



c 



I) 



+ 



, (15) 



where the external gluons carry soft momenta. Dots inside the brackets refer to terms with 
powers of 1/k different from 4 (i.e. corresponding to higher dimension operators or to identity 



operator n). From eq. (15) we get 



. (g^\A r {k)A s a (-k)\g b u ) 
(g a M 2 \g b ) 



'pa 



kpkcr 



Z 3 (M a )Z7 3 V)c2(^,«(M)) +••• 



k 2 



= Z- L (fi 2 ) c 2 ( j^,a(p) ) +... , (16) 



where Ar = Z 3 l ^ 2 A and A\ = Z~^A 2 , while Z = Z^ l Zj^2 is a useful notation denoting the 
divergent factor of the matrix element coming from proper vertex corrections. If one takes the 
logarithmic derivatives with respect to /x in both Eq. (|i~6|)'s hand-sides, the following differential 
equation is obtained: 



d d ) f k 2 

-2 7 (a(n)) A 2 + 2 7 (a(/i)) + + p (a(/i)) — j c 2 ( -j , a(/i) 







(17) 



An analogous differential equation describing the behaviour of the three-point Wilson coef- 
ficient on the renormalization momentum, fi, can be obtained similarly, 



d d } (k 2 

-2 7 (a(/x)) A 2 + 37 (a(/x)) + + P («(m)) ^ j d 2 ( -5 , «(m) 



(18) 



Here we have defined : 



7A2 



d In ft 2 



InZ 



A 2 



(19) 



3 It should be remembered that other terms, like d^A^, with the same dimension of A 2 , do not survive 



and 7 (a(fj,)) is the gluon propagator anomalous dimension. Reexpressing these evolution equa- 
tions in terms of 

7(«(M)) = ^^InZ(^) . (20) 

We obtain : 

5 Tn7 +/3(a( ^^} C2 U 

and 

The leading log solution for both Eqs. (21,22) can be written as, 



2 7 («(/*)) + thtt: + P («(/*)) — !• c 2 ( 1=0; (21) 



27 (aOi)) + 7 (a0*)) + T^TT: + («M) — \ d 2 [ ~ > «(/*) ) = ■ (22) 



^ 2 \A ( az(k)\~'Po 



c 2 -2,a(/x) = c 2 (l,a(A;)) . 

"fQ — 2"YQ 

d 2 Qj, a(Ai)) = cfe(l,a(fc)) (^|) ^ • (23) 

where 70 is defined in analogy with 70 : 

7(a( M )) = -7 ^gl + ... (24) 

The prefactors c 2 (1, a(/c)) and d 2 (1, ce(k)) have to be matched at tree level to Eq. (|TT|) . The 
only solutions are of the form : 



c 2 (!,«(&)) = 3 1 + 



1 



log(/c/A QC D) 



rf2(1 , aW)= 9^)(l + ( log(fc/ 1 AQcD) || . (25) 



The C ^ io g (fc/AQCD) J ^ erms are clearly of the same order as the next-todeading contributions 
to the anomalous dimension which are systematically omitted in this paper. 

Of course, these solutions of Eqs. (21,22) define the dependence of the Wilson coefficients not 



only on the renormalization momentum, /x, but simultaneously on the momentum scale k 2 . This 
is a straightforward consequence of standard dimensional arguments: the only dimensionless 
quantities are the ratio k 2 / fi 2 and a. Then, as soon as one knows perturbatively 7 («(//)), 
7(a(/x)) and (3(a(fi)), the leading logarithmic behaviour on k is available. 

As already mentioned, the gluon propagator anomalous dimension and the beta function are 
known up to three loops in the MOM scheme and up to three and four loops, respectively, in the 
MOM. The anomalous dimension of the A 2 operator is obviously less stimulating for perturba- 
tive QCD community. We have done this calculation to one loop (see appendix), obtaining 



and 



7 (a( M )) = -70^ + - = -—-4^- + - • (26) 
, , ,, d , „ 35A r C 7 01(a) 

™<°M>-3ta? lnZ "--na--r + " («) 
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2.3 Gluon propagator with leading logs for the condensate coefficient 

Let us now specify our approach to lattice results. Using the definitions in Eqs. (0) and (| 
we will match our lattice results to 



k 2 Gf(k 2 ,n 2 ) + Q(a 2 ) , (28) 



Latt 

where the adequate control of lattice artifacts^ reduces the UV discretization errors to an ac- 
ceptable level. From Eq. (P), 



2 



'k 



c 2 73,a(A») 



where we explicitly factorize the Wilson coefficient of the identity operator which, as was previ- 
ously indicated, is known to three loops. Nevertheless, for consistency, all the terms inside the 
parenthesis in the r.h.s of Eq.(]29|) will be developed only to leading order, including cq ( t-, «(//)), 



7Q 

k 1 \ / a(fc)\ 0o 



co ' Lo U ,a(/i) J = w^yJ (30) 

Terms of the order of 0(l/(k 2 /i 2 )) have been neglected, as well as, of course, those of 0(l/k 4 ) 
coming from higher dimension operators. One free parameter, i.e. a boundary condition, has 
to be fitted from lattice data. It can be either a(/i) or the A parameter, i.e. the position of 
the perturbative Landau pole. We choose the latter. We write cfynoop i n terms of the MOM 
coupling constant^] and the A parameter in Eq. ( |30|) in the MOM scheme^], 

A = A MOM ~ 3.334 . (31) 

We finally obtain 

Z^T(k 2 ,a) = Z L M rV,a)co(^,«(/i) 

TO +70 1 

where 



n 2 7Q +7Q 

R {2) = 5 Z T\ P ° (33) 



Po(Ni-l) V A 



4 See ref. where we discuss at length the artifacts of the lattice gluon propagator evaluation. 
5 MOM, for instance, or whatever renormalization scheme could be used alternatively. Our preference for the 
MOM scheme has been ex pla ined above. 
6 See, for instance, ref. H13I] 



2.4 Running coupling constant 



By taking the OPE expansions in Eqs. (10) and @ (0) can be written as 



9R(k 2 



with the identification 



5R, P ert(fe 2 ) < 1 



4(JV*-1) A; 2 




^2 



4 ( ^, 



(34) 



A; 2 

5R,pert(fe 2 ) = <Z (-2,a(//) 



-3/2 



(35) 



Applying the results given by Eqs. ( 23]) 32-35]), «mom = 5r,/(47t) verifies 



^ 70 +70 1 

'WlOM(fe) = apcrt(k) (l + -R (3) ( ln ^j ^2 



(36) 



where 



18vr 2 



00 (iV|-l) V~A 
Again, we do not retain 0(l/(k 2 fj, 2 ), 1/k )-terms 



7Q +7Q 



(37) 



3 Fitting the data to our ansatze 

We shall follow in this section the OPE testing approach proposed in ref. Q: trying a 
consistent description of lattice data for two- and three-gluon Green functions from ref. f|,[Zl, 
||. We are however in a much better position than in H. In the latter work, only two-loop 
information was available for the beta function and do while the subdominant Wilson coefficients 
C2 and d.2 were computed only at tree-level. This had the practical inconvenience of preventing 
a simultaneous fit of both A MS and (A 2 ): the A MS parameter had to be taken from outside our 
matching procedure. Now the new input for three-loop MOM beta function and Co coefficient Q 
enable us to perform a self-consistent test by combining the matching of the gluon propagator 
and of «MOM 



to formulas in Eqs. (32- p3jj36H37|) , where the three quantities, Aj^g- and gluon 
condensates from both Green functions, are taken to be fitted on the same footing. Of course, 
the test consists in checking the equality of the two gluon condensates obtained from those two 
different Green functions. 

In the case of the gluon propagator, the factor co(l, a(/i))Z^^ t M (/i 2 , a), which carries all 
the logarithmic dependence on the lattice spacing, appears as an additional parameter to be 



fitted. As we explained in ref. [0], a large fitting window is an important "ace" to restrict the 
potentially dangerous confusion between Wilson coefficients for different powers. To combine 
data over such a large energy window we need to match the lattice results obtained with different 
lattice spacings and the last factor carrying lattice spacing dependence should be independently 
fitted for each one. On the contrary, the running coupling constant should be regularization 
independent and the matching of data sets corresponding to different lattice spacing can be 
imposed without tuning any additional parameter (this is by itself a positive test of the goodness 
of the procedure used to build our data set). As a matter of fact, this is why the matching 
of the latter to perturbative formulas is much more constraining than that of the former in 
order to estimate Aj^g, as discussed in refs. || (/], H) • The details of the lattice simulations, of 
the procedures used to obtain an artifact-safe data set or of the definition of regularization- 
independent objects permitting lattice regularized data to be matched to continuum quantities 
in any scheme, can be found in those references. We will now present the results of the fitting 
strategy just described. 



Two loop fit We first perform the combined fit for two- ant three-gluon Green function at 
two-loop level for the leading Wilson coefficients. In fig. |2], we plot lattice data and the curves 
given by Eqs. (32-^3,3^-37) with the following best-fit parameter: 



propagator : J (A 2 ) R ^ = 1.64(17)GeV; a 



mom . 



(A^)r^ = 3.1(3)GeV. 



alpha 



1.86(4); = 172(15)MeV 



prop 



(38) 



with a x 2 /d-o.f ~ 1.1 for the combined fit. 
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Figure 2: Comparison of the 2-loops fit to the ratio of the renormalization constants at k and at 10 GeV 
and to a s (k) with the lattice data for 2.5 < k < 10 GeV . The dotted line shows the perturbative part. 
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Three loops fit The present perturbative knowledge allows a three-loop level fit for leading 
Wilson coefficients. Analogously to the previous paragraph, we plot in Fig. |3] the lattice data 
and curves given by Eqs. (|32|- 33|,36-37), -^pe^ an d a pe rt taken at three loops, with the following 
best-fit parameters: 



propagator : 



alpha 



(A 2 ) R ^ = 1.55(17)GeV; Q • ^(A 2 ) R ^ = 1.9(3)GeV. 

1.21(18); A m = 233(28)MeV (39) 



prop 



with x 2 /d.o.f ~ 1.2. Combining the results obtained from a M M and from the propagator in 
the standard way gives our final result W (A 2 ) R a = 1.64(15)GeV. The renormalization scale [i is 
taken to be 10 GeV in both combined fits at two- and three-loop level. However we have checked 
that, varying ^ over the fitting window where we can legitimately neglect terms of 0(l/(/j, 2 k 2 )) 
in Eq. (|37l), the ratios of condensates in Eqs. (3J,3^) remain essentially unmodified. In fact, 
that in Eq. (32) does not depend on [i has been explicitly tested over the fitting window 
(the same is obvious for in Eq. (1361) where nothing depends on fj). 
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Figure 3: Same as Figure 2 at 3-loops level. 



4 Discussion and conclusions 

The gluon condensate (A 2 ) has been computed from the deviations of both the lattice non- 
perturbative evaluation of the MOM as and the gluon propagator from their known perturbative 
behaviour. We have described these non-perturbative deviations using OPE, and fitting the 
condensate to match both sides. The use of the self-consistent fitting strategy described in the 
preceding section leads simultaneously to a prediction A^jg and to two independent estimates 
of(^l 2 ). 

The fit using the two-loops perturbative expressions for both the MOM as and the gluon 
propagator clearly fails: a clear disagreement between the two independent estimates of w (A 2 ) 

1 1 



is found. The ratio of both estimates is 1.86(4) from Eqs. (|38D . This confirms the preliminary 
analysis in ref. Q, where only tree-level Wilson coefficients were computed. In this preliminary 
work, a self-consistent three-loop analysis was not possible because the MOM beta function 
was not known up to three-loops. Nevertheless we tried to fit the third coefficient of the beta 
function, to reach a good agreement between the two estimates of (A 2 ), the Aj^g parameter 
being taken from previous works to be the same for both two- and tree-point Green function 
matchings. Our estimate, (3% = 7400(1500) was about twice larger than the result [3% = 3088 in 

0. Still this fit went in the right direction, whence the authors of Q expected their result to 
lead to a fair fit to lattice data. 

This expectation turns out to be correct : 

First, the ratio of the two estimates of W (A 2 ) is equal to 1.21(18), i.e. compatible with 

1, provided the leading Wilson coefficients are consistently expanded at the three-loop level and 
the subleading coefficients of (A 2 ) are computed to the leading logarithms. Second, in the same 
joint fit, Aj^jg is estimated to be 233(28)Mey, in amazing agreement with previous estimates 
of Ajyfg- appearing in the literature (see for instance Thus, the present analysis ends 
up with a twofold success and we can conclude that OPE leads to a good description of the 
deviations of the running coupling constant and of the gluon propagator from their perturba- 
tive behaviour in terms of perturbatively available coefficients multiplying one phenomelogical 
condensate: the sole non vanishing non-perturbative contribution up to the order 1/p 2 , namely 
the gluon condensate (A 2 ). 

However, for this OPE description to be consistent, it is unambiguously demonstrated that 
the leading coefficient must be taken to three loops; on the contrary there is a clear failure at 
two loops, even though our analysis has been performed at an unusually large energy scale, up 
to 10 GeV. As was discussed in the preliminary study, such a disagreement indicates that we 



are, in this case, in the situation described in [11|, i.e. that the perturbative order is too low to 



give an acceptable precision in the estimate of the power corrections. Taking into account power 
corrections does not make sense unless the leading contribution is computed perturbativeley to 
a sufficient accuracy. In a simple mathematical language, it makes no sense to consider the 1 jp 2 
corrections when one does not consider the l/ln(p 2 ) corrections to a sufficient order. OPE is 
often used with the leading coefficients only known to two loops (sometimes to one loop). We 
believe in view of our results that these attempts should be reconsidered with care. 

As for us, we were in a particularly favourable situation to analyse the problem thoroughly. 
We have rather accurate results. The dimension 2 power correction clearly shows up clearly 
and can be fully consistently attributed to an A^A^ condensate in the Landau gauge in full 
agreement with the theoretical expectations . Since we are working in the Landau gauge we 
produce (and use later on) bare gauge field configurations which minimize A^A^ with respect 
to the gauge group. This has the interesting consequence that the quantities we measure are 
invariant under infinitesimal gauge transformations in the vicinity of the Landau gauge. Still 
the link between what we call the (A 2 ) and the (^^in) defined in || should be better clarified, 
which implies a better understanding of the renormalization proceduref]. Taking such a direct 
link for granted, we can estimate from Eqs. (^) the tachyonic gluon mass defined in ref. Q, to 
be ~ .8 GeV^. Using instead the notion of critical mass, M 2 lit , introduced in ref. |2C], which is 



the scale at which the non-perturbative condensate contributes 10% of the total, we estimate it 
for the gluon propagator to be ~ 2.6 GeV. Both these scales express a rather large contribution 
from the A 2 condensate. 



7 it is a pleasure to acknowledge discussions with V.I Zakharov on this topic 
8 We recall that all scale-dependent quantities are evaluated at fi — 10 GeV. 
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A One-loop anomalous dimension of A 2 

The task of computing the one- loop anomalous dimension of the matrix element (A 2 ) (that 
of the local operator itself is directly obtained from the former as explained in section § ) 
requires only to isolate the UV divergent part of the diagrams in Fig. ^. We follow dimensional 
regularization prescriptions to write: 

« (9, -1) = (sT ~ ^f) ^ Q {3^2 + O K 2 ) }+•••) , 

ME -1) = (f - S<* Q {-lNo% + O (of) I +...); (40) 

where an a(b)-like diagram in Fig. |2| with amputated external gluon legs is denoted by ^ a (b) an d 
where e = 2 — d/2. The particular kinematics we choose (see the figure), where the incoming 
momentum flow is non vanishing, eliminates automatically the IR divergences and makes the 
UV analysis easier. We require the final result for amputated diagrams to be transverse to the 
external momenta, but this is merely a convention to be also applied to the tree-level term. 
Furthermore, had we considered two different incoming momenta, the tensors in Eq. (|40|) would 
have acquired a more complicated form. 




q+k 



Figure 4: The graphs involved in the computation of the anomalous dimension of A 2 . The cross-hatched 
blobs indicate the insertion of the A 2 operator, the dots are ordinary QCD vertices. 

This kind of IR regularization, by imposing a non-null incoming momentum flow to the local 
operator, leads of course to UV poles results equivalent to those obtained from any other one. We 
have tested this by considering a null incoming momentum flow and both introducing a certain 



cut-off to regularize IR divergences and separating IR and UV poles obtained by dimensional 
regularization. 

If we collect the tree- level results and those from Eqs. (|40|), we can write: 



«HV> = 2 - S*) J- (i + i + (<*)} + ■ ■ ■) < 41 > 

where the matrix element in l.h.s. of Eq. ( |4l| ) is defined for explicitly cut external gluons. 
Combinatorics gives a multiplicity factor 2 for a-like diagram, 1 for b-like, which have been 
taken into account in the last result. 

We should now renormalize the matrix element defined in Eq. (|4l]). Our aim being to 
determine its anomalous dimension computation to only one-loop, the simple MS prescription 
of simply dropping away from bare quantities the poles for e — > can be applied. Discrepancies 
between such a prescription and MOM or any other one appear only beyond one-loop. Then we 
will have : 



(42) 



From Eq. (42), the anomalous dimension can be written as (see for instance ||12| ): 



7 MS (* MS M) 



din fi 2 



InZ 



MS 



3N C a MS {fi) 

4 47T 

Thus we obtain in the MOM scheme, 



(43) 



7(«<tf) = -^^ + 0( ( a(^)) 2 
From this we deduce finally, including the gluon anomalous dimension 

35N C a(ji) 



j A 2 (a(p)) 



12 



4tt 



+ Ol(a(jj))' 
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